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Abstract 

We describe the volume dependence of matrix elements of local fields to all orders in 
inverse powers of the volume (i.e. only neglecting contributions that decay exponentially 
with volume). Using the scaling Lee- Yang model and the Ising model in a magnetic field 
as testing ground, we compare them to matrix elements extracted in finite volume using 
truncated conformal space approach to exact form factors obtained using the bootstrap 
method. We obtain soHd confirmation for the form factor bootstrap, which is different 
from all previously available tests in that it is a non-perturbative and direct comparison 
of exact form factors to multi-particle matrix elements of local operators, computed from 
the Hamiltonian formulation of the quantum field theory. We also demonstrate that 
combining form factor bootstrap and truncated conformal space is an effective method 
for evaluating finite volume form factors in integrable field theories over the whole range 
in volume. 



1 Introduction 

The matrix elements of local operators (form factors) are central objects in quantum held 
theory. In two-dimensional integrable quantum held theory the S matrix can be obtained ex- 
actly in the framework of factorized scattering (see [1, 2] for reviews), and using the scattering 
amplitudes as input it is possible to obtain a set of axioms satished by the form factors [3], 
which provides the basis for the so-called form factor bootstrap (see [4] for a review). 

Although the connection with the Lagrangian formulation of quantum held theory is rather 
indirect in the bootstrap approach, it is thought that the general solution of the form factor 
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axioms determines the complete local operator algebra of the theory [5]. This expectation was 
confirmed in many cases by explicit comparison of the space of solutions to the spectrum of 
local operators as described by the ultraviolet limiting conformal field theory [6, 7, 8, 9]; math- 
ematical foundation is provided by the local commutativity theorem stating that operators 
specified by solutions of the form factor bootstrap are mutually local [4]. Another important 
piece of information comes from correlation functions: using form factors, a spectral repre- 
sentation for the correlation functions can be built which provides a large distance expansion 
[10, 11], while the Lagrangian or perturbed conformal field theory formulation allows one to 
obtain a short-distance expansion, which can then be compared provided there is an overlap 
between their regimes of validity [11]. Other evidence for the correspondence between the 
field theory and the solutions of the form factor bootstrap results from evaluating sum rules 
like Zamolodchikov's c-theorem [12, 13] or the A-theorem [14], both of which can be used to 
express conformal data as spectral sums in terms of form factors. Direct comparisons with 
multi-particle matrix elements are not so readily available, except for perturbative or 1/N 
calculations in some simple cases [3]. 

Therefore, part of the motivation of this paper is to provide non-perturbative evaluation 
of form factors from the Hamiltonian formulation, which then allows for a direct comparison 
with solutions of the form factor axioms. Another goal is to have a better understanding of 
finite size effects in the case of matrix elements of local operators, and to contribute to the 
investigation of finite volume [15] (and also finite temperature [16]) form factors and correlation 
functions. 

Based on what we learned from our previous investigation of decay rates in finite volume 
[17], in this paper we determine form factors using a formulation of quantum field theory in 
finite volume. In two space-time dimensions this is most efficiently done using the truncated 
conformal space approach (TCSA) developed by Yurov and Zamolodchikov [18], but one could 
have also made use of, for example, lattice field theory: the only important point is to have a 
method which can be used to determine energy levels and matrix elements of local operators as 
functions of the volume. We give a relation between finite and infinite volume multi-particle 
form factors, which is a natural extension of the results by Lellouch and Liischer for two- 
particle decay matrix elements [19]. It is used to determine form factors in two important 
examples of integrable quantum field theory: the scaling Lee- Yang model and the Ising model 
in a magnetic field. We show that these results agree very well with the predictions of the form 
factor bootstrap. We only treat matrix elements of local fields between multi-particle states 
for which there are no disconnected pieces (which appear whenever there are particles with 
coincident rapidities in the left and right multi-particle states); the treatment of disconnected 
pieces, together with a number of theoretical arguments and ramifications, are postponed to 
a subsequent publication [20]. 

The organization of the paper is as follows. In section 2, after a brief review of necessary 
facts concerning the form factor bootstrap we give the general description of form factors 
(without disconnected pieces) to all orders in 1/L based on an analysis of two-point correlation 
functions. Section 3 describes the two models we chose for demonstration, and also specifies 
the method of evaluating form factors from truncated conformal space, Numerical results on 
elementary form factors (vacuum-many-particle matrix elements) are given in section 4, while 
the general case is treated in section 5. We give our conclusions in section 6. 
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2 Form factors in finite volume 



2.1 Form factor bootstrap 

Here we give a very brief summary of the axioms of the form factor bootstrap, in order to set 
up notations and to provide background for later arguments; the interested reader is referred 
to Smirnov's review [4] for more details. Let us suppose for simplicity that the spectrum of the 
model consists of particles Ai, i = I, . . . , N with masses rrii, which are assumed to be strictly 
non-degenerate, i.e. rrii 7^ rrij for any i 7^ j (and therefore also self-conjugate). Because of 
integr ability, multi-particle scattering amplitudes factorize into the product of pairwise two- 
particle scatterings, which are purely elastic (in other words: diagonal). This means that 
any two-particle scattering amplitude is a pure phase, which we denote by Sij (6) where 6 is 
the relative rapidity of the incoming particles and Aj. Incoming and outgoing asymptotic 
states can be distinguished by ordering of the rapidities: 



1^1, 



and states which only differ in the order of rapidities are related by 
1^1, ...,9k, Ok+l, ■ ■ ■ , dn)ii...ikik+i-in = Si^,if.j^^{6k — ^fe+l)|^l5 • • • > dk+l,9k, ■ ■ ■ , 9n)ii...ik+iik--in 

from which the S matrix of any multi-particle scattering process can be obtained. The nor- 
malization of these states is specified by the following inner product for the one-particle states: 

j{e'\e)i = Sij27rS{9' -6) 

The form factors of a local operator 0(t,x) are defined as 

-^mn(^l> • • • >^ml^l> • • • > ^n)ji...j^;n...i„ = ji...jm{^l, ■ ■ 6'm|C(0, 0)1611, • • ■,dn)h...in (2-1) 

With the help of the crossing relations 

-^mn(^l5 • • • ' ^ml^l' ■ ■ ■ ' ^n)ji...jm;h—in ~ 

^m-\n+l{^lt ■■■ 1 ^m-l\^m + ^TT, ^l, . . . , On)ji...jm-i;jmil---in 

n k—1 



k=l 1=1 
-^m-ln-l(^l) • • • ) ^m-ll^l5 • • • ) ^fc-li (^k+1 ■ ■ ■ , 0n)ji...jm~i;jmil-ik-lik+l---'>'n^ (2-2) 

all form factors can be expressed in terms of the elementary form factors 



F^{9i, 9n)n...ir. = (0P(0, 0)1^1, ... , 9n) 



11. ..In 



which satisfy the following axioms: 
I. Exchange: 



{9i, . . . , 9k, 9k-{-l, • • • , 9n)ii...ikik+i---in ~ 
Sikik+ii^k - 0k+l)F^{9i, . . .,9k+l,9k, . . . ,9n)h...ik+iik-in (2-3) 
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II. Cyclic permutation: 

F^{ei + 2m,e2,...,en) = FS{e2,...,en,ei) (2.4) 

III. Kinematical singularity 



-i Res Ff+2(^ + ^i, • • • , ^n)i,n...i„ =1-5^^-11 " ^'^) ^"""(^i' • • • ' 

V fe=i / 



(2.5) 



IV. Dynamical singularity 

-i Res F„0+2(^ + e' - ivij2, Oi,..., en)ijn...in = ^i, • • • , On)kH...i„ (2.6) 

9=9 

whenever k occurs as the bound state of the particles i and j, corresponding to a bound state 
pole of the S matrix of the form 

S,j{e ^ iu'ij) ^ -^^^ (2.7) 

where F^^- is the on-shell three-particle coupling and u^j is the so-called fusion angle. The 
fusion angles satisfy 

'^fc = mf + m'j + 2mimj cos u^j 

and we also used the notation ufj = tt — n^^-. Axioms I-IV are supplemented by the assumption 
of maximum analyticity (i.e. that the form factors are meromorphic functions which only 
have the singularities prescribed by the axioms) and possible further conditions expressing 
properties of the particular operator whose form factors are sought. 

2.2 Finite size corrections for form factors 

Let us consider the spectral representation of the Euclidean two-point function 
(O(x)O'(0,0)) = E E n / ^ if (^i,02,...,0n)n..... X 

if (01,02,..., ^n)+...i„cxp ( -rEmi,cosh0fc ) (2.8) 



k=l 



where 



F^' {01,02, 0n)Z...i^ = 0n\O'iO, 0)|0) = i^f (01 + in, 0^ + in, . . . ,0^ + in)i,..,„ 

(which is just the complex conjugate of F^ for unitary theories) and r = Vt^ + is the 
length of the Euclidean separation vector x = (r, x) . 
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In finite volume L, the space of states can still be labeled by multi-particle states but 
the momenta (and therefore the rapidities) are quantized. Denoting the quantum numbers 
Ii, . . . ,In the two-point function of the same local operator can be written as 

oo 

(O(r,0)O'(0,0))z, = EE E ^0|^(0'0)K^i'^2,...,/n})n....„,LX 

n=0 ii...i„ h...In 

h...ij{h,l2, . . . , In}\O'{0, 0)|0)l exp ^-r ^ rm, cosh^fe^(2.9) 

where we supposed that the finite volume multi-particle states h, ■ ■ ■ , In})ii...in,L are 
orthonormal and for simplicity restricted the formula to separation in Euclidean time r only. 
The index L signals that the matrix element is evaluated in finite volume L. Using the finite 
volume expansion developed by Liischer in [21] one can easily see that 

(0(r, 0)0' (0)) - (0(t, 0)O'(0))l ~ 0(e-'^-^) (2.10) 

where /x is some characteristic mass scale. 

Note that Liischer's finite volume expansion is derived using covariant perturbation theory, 
but the final expansion is obtained after resummation in the coupling constant, and is expected 
to hold non-perturbatively. The finite volume corrections are then expressed in terms of the 
exact one-particle-irreducible vertex functions of the infinite volume theory. According to 
Liischer's classification of finite volume Feynman graphs, the difference between the finite 
and infinite volume correlation function is given by contributions from graphs of nontrivial 
gauge class, i.e. graphs in which some propagator has a nonzero winding number around the 
cylinder. Such graphs always carry an exponential suppression factor in L, whose exponent can 
be determined by analyzing the singularities of the propagators and vertex functions entering 
the expressions. In a massive theory, all such singularities lie away from the real axis of the 
Mandelstam variables, and the one with the smallest imaginary part determines /i. It turns 
out that the value of is determined by the exact mass spectrum of the particles and also 
the bound state fusions between them [21, 22]. Therefore it is universal, which means that 
it is independent of the correlation function considered. In general fj, < mi where mi is the 
lightest particle mass (the mass gap of the theory), because there are always corrections in 
which the lightest particle loops around the finite volume L, and so the mass shell pole of the 
corresponding exact propagator is always present. Contributions from such particle loops to 
the vacuum expectation value are evaluated in subsection 4.1 (for a graphical representation 
see figure 4.2), while an example of a finite volume correction corresponding to a bound state 
fusion (a so-called /u-term) is discussed in subsection 4.1.2 (figure 4.7). 

To relate the finite and infinite volume form factors a further step is necessary, because 
the integrals in the spectral representation (2.8) must also be discretized. Let us consider this 
problem first for the case of free particles: 




where 

Pk = "life sinh Ok , ujk = ruik cosh 9k 
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are the momenta and energies of the particles. In finite volume 

27rJjfc 

and it is well-known (as a consequence of the Poisson summation formula, cf. [23]) that 

5 • • • ' ^) = (^) " (n dp,^ siPi, . . . ,Pn) + 0{L-^) (2.11) 

provided the function g and its first N derivatives are integrable. Recalling that form factors 
are analytic functions for real momenta, in our case this is true for derivatives of any order, 
due to the exponential suppression factor in the spectral integrals, provided the form factors 
grow at most polynomially in the momentum, i.e. 

|i^n(6'i + ^,^2 + ^,...,^n + 6')| -e^l^l as |6'|^oo 

This is true if we only consider operators which have a power-like short distance singularity 
in their two-point functions [24]: 

(0|O(x)O(0)|0) = ^ 

Such operators are called scaling fields and are generally assumed to form a closed algebra 
under the operator product expansion. 

Therefore the discrete sum differs from the continuum integral only by terms decaying 
faster than any power in 1/L, i.e. by terms exponentially suppressed in L. Taking into 
account that (2.8) is valid for any pair O, O' of scaling fields, we obtain 

(0|O(0, 0)|{/i, . . . , /n})n...^„,L = , ^ =FS{h, 9n)n...ir. + ©(e"'^'^) (2.12) 



where 



and 



sinh 9k 



n 



{61,..., 9n) = Ylmi^L cosh Ok (2.13) 



k=l 



Pn^ is nothing else than the Jacobi determinant corresponding to changing from the variables 
27rlfc to the rapidities dk- The term 0{e~^ ^) signifies that our considerations are valid to all 
orders in 1/L, although our argument does not tell us the value of 11': to do that, we would 
need more information about the correction term in the discretization (2.11). 

In the case of interacting particles a more careful analysis is necessary because the quan- 
tization rules are different from the free case. In a two-dimensional integrable quantum field 
theory, general multi-particle levels are determined by the Bethe-Yang equations 

(5fe(6*i,...,^n)n...i„ ="^jfc^sinh6ife + E'^«fe'((^fc~^') = , k = l,...,n (2.14) 

l^k 
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which are valid to any order in 1/ij, where 



5iji9) = -ilogS^j{e) (2.15) 
are the two-particle scattering phase-shifts. Therefore the proper generalization of (2.12) is 

(0|O(0, 0)|{/i, . . . , In})^,...^„,L = [ =F^i9i, 9n)n...ir. + ©(e"'^'^) (2.16) 

where 



Pn...iAGi,...,0n) = detj^''\ei,...,dn)n...ir. (2.17) 

r{n)fp a ^ _ dQkjdi,. . . ,9n)h...ir, 

39, 



^kl V'iT ■ ■ ^^n)i-i...in — ) /C, t — 1, . . . , fi 



and 9k are the solutions of the Bethe-Yang equations (2.14) corresponding to the state with 
the specified quantum numbers /i, . . . , I„ at the given volume L. 

Similar arguments were previously used to obtain the finite size dependence of kaon decay 
matrix elements by Lin et al. [25]; they also give a more detailed analysis of the discretization 
in (2.11). In principle it could be that /i' < /v, if the discretization introduces larger errors 
than the finite size correction of the two-point functions. However, it is possible to give an 
argument that actually fi' = (i. Recall that the Poisson formula gives the discrete sum in terms 
of a Fourier transform: the leading term is the Fourier transform of the summand evaluated 
at wa\'e number (i.e. the integral) and the corrections are determined by the decay of the 
Fourier transform at large wave numbers. The function we need to consider is 



i ^ 



ii...in 

exp ^-r rriii^ cosh 9kj (2.18) 

where Pk = mii^smh.9k are the momentum variables. Due to the analyticity of the form 
factors for real rapidities, this function is analytic for physical (real) momenta, and together 
with all of its derivatives decays more rapidly than any power at infinity. Therefore its Fourier 
transform taken in the momentum variables has the same asymptotic property, i.e. it (and its 
derivatives) decay more rapidly than any power: 



h{Ki, ...,«;„) ~ e 



for large k. As a result, discretization introduces an error of order c^^"^. The asymptotic 
exponent /x" of the Fourier transform can be generally determined by shifting the contour 
of the integral transform and is given by the position of the singularity closest to the real 
momentum domain (this is essentially the procedure that Liischer uses in [21]). Singularities 
of the form factors are given by the same analytic structure as that of the amplitudes which 
determine the exponent fi in eqn. (2.10) (see the discussion there). Thus we expect that 
IJ," = 11 and, as a result, n' = ji. 



7 



This argument is just an intuitive reasoning, although it can be made a httle more precise. 
First of all, we must examine whether the determinant Pij^...i„{Oi, . . . ,9n) can have any zeros. 
It can always be written in the form 

Pn...iA<^i,...,0n)= (^f[m,Lcosh9k^ (1 + 0(1/L)) 

The leading factor can only be zero when 0^ = ^ for some k, which corresponds to Pk = irrik, 
giving fx" = rrik in case this is the closest singularity. That gives a correction 



which is the same as the contribution by an on-shell propagator wound around the finite vol- 
ume, and such corrections are already included in the e"'^^ term of (2.10). Another possibility 
is that some phase-shift function 6(6) in the 0{1/L) terms contributes a large term, which 
balances the 1/L prefactor. For that its argument must be close to a singularity, and then 
according to eqn. (2.7) we can write 

6{0) ~ log(^ - Tx) ~ 0{L) 

where u is the position of the singularity in the phase-shift^. This requires that the singularity 
is approached exponentially close (as a function of the volume L), but the positions of all 
these singularities are again determined by singularities of the vertex functions, so this gives 
no new possibilities for the exponent /v/'. 

A further issue that can be easily checked is whether the Fourier integral is convergent for 
large momenta; the function (2.18) is cut off at the infinities by the factor 

n 

JJ exp(-mij.r cosh Ok) 
k=l 

which can only go wrong if for some k 

SRe cosh 6k <0 

but that requires 

9m6'fe > - 

which is already farther from the real momentum domain then the position of the on-shell 
propagator singularity. 

Obtaining a sound mathematical proof requires applying Liischer's finite volume expansion 
directly to form factors which is out of the scope of the present work. The corrections computed 
in subsection 4.1 can be considered as the simplest examples, but we intend to take this line 
of investigation further in the future [26] . 

^Strictly speaking, this argument only works in an integrable field theory where we know the form (2.17) 
of the determinant, and also know the singularity structure of the phase-shift S: the two-particle S matrix 
amplitude as a function of the rapidity can only have poles (possibly of higher order, but that does not affect 
the main line of reasoning here). 
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We also remark that there are no finite volume states for which the quantum numbers of 
any two of the particles are identical. The reason is that 

SiiiO) = -1 

(with the exception of free bosonic theories) and so the wave function corresponding to the 
appropriate solution of the Bethe-Yang equations (2.14) vanishes. We can express this in 
terms of form factors as follows: 



(0|O(0,0)|{/i,/2,...,/n})*i...i„,L = 

whenever Ik = Ii and ik = h for some k and /. Using this convention we can assume that the 
summation in (2.9) runs over all possible vahies of the quantum numbers without exclusions. 
Note that even in this case the relation (2.16) can be maintained since due to the exchange 
axiom (2.3) 

Pni^l^ ■ ■ ■ )^n)n...i„ = 

whenever 9k = 9i and ik = ii for some k and I. 

It is also worthwhile to mention that there is no preferred way to order the rapidities on the 
circle, since there are no genuine asymptotic in/out particle configurations. This means that 
in relation (2.16) there is no preferred way to order the rapidities inside the infinite volume 
form factor function Ff^ . Different orderings are related by ^'-matrix factors according to the 
exchange axiom (2.3), which are indeed phases. Such phases do not contribute to correlation 
functions (cf. the spectral representation (2.8)), nor to any physically meaningful quantity 
derived from them. In subsection 4.2.1 we show that relations like (2.16) must always be 
understood to hold only up to physically irrelevant phase factors. 

The quantity (^i, • • • , On) is nothing else than the density of states in rapidity space. 

It is also worthwhile to mention that relation (2.16) can be interpreted as an expression for 
the finite volume multi-particle state in terms of the corresponding infinite volume state as 
follows ^ 

\{Il,...,In})ii...i„,L = I _ _ \dl,---,dn)ii...in (2-19) 

Y Pn...jn(^l) ■■■,Gn) 

This relation between the density and the normalization of states is a straightforward appli- 
cation of the ideas put forward by Saleur in [27]. Using the crossing formula (2.2), eqn. (2.19) 
allows us to construct the general form factor functions (2.1) in finite volume as follows: 

l'm}\OiO, 0)|{/i, . . . , /„})n...«„,L = 
^m+ni^'m + iTT, ■■■ ,0'i + 17^,01, ... , On) j^...j^n...i„ _^ Qf^^-^iL-j ^3 20) 

\/ Ph...in{Ol, On)Ph-jm{K-> ■■■■>^'m) 

provided that there are no rapidities that are common between the left and the right states 
i.e. the sets |6'i, . . . and . . . are disjoint. The latter condition is necessary to 

eliminate disconnected pieces. 

We stress that eqns. (2.16, 2.20) are exact to all orders of powers in 1/L; we refer to the 
corrections non-analytic in 1/L (eventually decaying exponentially as indicated) as residual 
finite size effects, following the terminology introduced in [17]. 
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3 Form factors from truncated conformal space 
3.1 Scaling Lee- Yang model 

3.1.1 Truncated conformal space approach for scaling Lee- Yang model 

We use the truncated conformal space approach (TCSA) developed by Yurov and Zamolod- 
chikov in [18]. The ultraviolet conformal field theory has central charge c = —22/5 and a 
unique nontrivial primary field $ with scaling weights A = A = —1/5. The cylinder of 
circumference L can be mapped unto the complex plane using 

27r , . , _ 2tt , /„ 1 N 

z = exp—[T — ix) , z = exp—[T + tx) (3.1) 
ij Ij 

The field $ is normalized so that it has the following operator product expansion: 

^{z, z)$(0, 0) = C(zz)^/''$(0, 0) + {zzf'H + ... (3.2) 
where I is the identity operator and the only nontrivial structure constant is 

C = 1.911312699- •• X i 
The Hilbert space of the conformal model is given by 

nLY= Vh®Vh 

h=0,-l/5 

where Vh (V/i) denotes the irreducible representation of the left (right) Virasoro algebra with 
highest weight h. 

The Hamiltonian of scaling Lee- Yang model takes the following form in the perturbed 
conformal field theory framework: 

jjSLY ^ jjLY ^-^ f ^x^Q^ ^) (3 3) 

Jo 

where ^ 

is the conformal Hamiltonian. When A > the theory above has a single particle in its 
spectrum with mass m that can be related to the coupling constant as [28] 

A = 0.09704845636 • • • x m^^/^ (3.4) 

and the bulk energy density is given by 

B = -§m' (3.5) 

The 5-matrix reads [29] 

^ sinh^ + zsin^ 

sinh 9 — isiiL^ 
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and the particle occurs as a bound state of itself at 9 = 27rz/3 with the three-particle coupling 
given by 

= -2^/3 

where the negative sign is due to the nonunitarity of the model. In this model we define the 
phase-shift via the relation 

Sly{0) = -e'^^^^ (3.7) 

so that (5(0) = 0. This means a redefinition of Bethe quantum numbers Ik in the Bethe-Yang 
equations (2.17) such they become half-integers for states composed of an even number of 
particles; it also means that in the large volume limit, particle momenta become 

msmhPfc = — - — 
ij 

Due to translational invariance of the Hamiltonian (3.3), the conformal Hilbert space H can 
be split into sectors characterized by the eigenvalues of the total spatial momentum 

P=^(Lo-Lo) 

the operator Lq — Lq generates Lorentz transformations and its eigenvalue is called Lorentz 
spin. For a numerical evaluation of the spectrum, the Hilbert space is truncated by imposing 
a cut in the conformal energy. The truncated conformal space corresponding to a given 
truncation and fixed value s of the Lorentz spin reads 

WTCs(s,ecut) = {iV') e W| {Lo-Lo)\^) =s\^), (Lo + Zo - ^) IV') = e|V) :e<ecut} 
On this subspace, the dimensionless Hamiltonian matrix can be written as 

ft, = ^(z,. + Io-^ + .(^G<-MB<-.) (3.8) 

where energy is measured in units of the particle mass m, I = niL is the dimensionless volume 
parameter, 

Gg^ = m (3.9) 
is the conformal inner product matrix and 

4'^= (^l^(^,^)|j)l.=.-=i (3.10) 

is the matrix element of the operator $ at the point z = z = I on the complex plane between 
vectors \j) from Hxcst'S, Ccut)- The natural basis provided by the action of Virasoro 
generators is not orthonormal and therefore G*-*-*^^ must be inserted to transform the left 
vectors to the dual basis. The Hilbert space and the matrix elements are constructed using 
an algorithm developed by Kausch et al. and first used in [30]. 

Diagonalizing the matrix hij we obtain the energy levels as functions of the volume, with 
energy and length measured in units of m. The maximum value of the cutoff Ccut we used was 
30, in which case the Hilbert space contains around one thousand vectors, slightly depending 
on the spin. 
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3.1.2 Exact form factors of the primary field $ 

Form factors of the trace of the stress-energy tensor were computed by Al.B. Zamolodchikov 
in [11], and using the relation 

e = zA7r(l-A)$ (3.11) 
we can rewrite them in terms of They have the form 



n n 



Fr,{er,...,en) = {^)HnQn{xi,...,Xn)\[ \{ l^^l-M. (3.12) 

i=l j=i+l ' J 

with the notations 

posn n — 1 



~ sinh f sinh ^ sinh ^ .^^ 



v{e) = exp f 2 / dt ^/r" ^ 7"" ^ e 
\ Jo tsmh TTt 



31/4 



\^2V2^;(0) J 

The exact vacuum expectation value of the field $ is 

($) = 1.239394325- • • x zm'^/^ 

which can be readily obtained using (3.4, 3.11) and also the known vacuum expectation value 
of e [11] 



(©) 



4^/3 

The functions Qn are symmetric polynomials in the variables Xi. Defining the elementary 
symmetric polynomials of n variables by the relations 

n n 

Yl{x + Xi) = J2 Vf ^ (xi , . . . , x„) , ) = for z > n 

i=l i=0 

they can be constructed as 

n ^ ^ (2) ^ (3) (3) 

Pn = det where M'lf = 4i-2j+i > i, J = 1, • • • , n - 3 

Note that the one-particle form factor is independent of the rapidity: 

= 1.0376434349 • • • x im.-^/^ (3.13) 
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3.2 Ising model with magnetic perturbation 

The critical Ising model is the described by the conformal field theory with c = 1/2 and has 
two nontrivial primary fields: the spin operator a with A^- = A^- = 1/16 and the energy 
density e with Ag = Ag = 1/2. The magnetic perturbation 



H 



= H^ + h [ dxaiO, 
Jo 



is massive (and its physics does not depend on the sign of the external magnetic field h). The 
spectrum and the exact S matrix is described by the famous Eg factorized scattering theory 
[31], which contains eight particles Ai, i = 1, . . . ,8 with mass ratios given by 

TT 

1712 = 2mi cos — 
5 

TT 

ma = 2mi cos — 
77r 

777-4 = 27772 COS — 
OU 

27r 

7775 = 27772 COS — 

15 

7776 = 27772 COS 

30 

TT 

7777 = 27774 cos — 

5 

TT 

7778 = 27775 COS — 

5 

and the mass gap relation is [32] 

mi = (4.40490857...) 1^1^/^^ 



or 



h = KhUif'^ , Kh = 0.06203236 . . . (3.14) 



The bulk energy density is given by 

B = -0.06172858982 ■ ■ ■ x mj (3.15) 
We also quote the scattering phase shift of two Ai particles: 

^/^x fll fll f2l ri sinh ^ + 7 sin TTX 

All other amplitudes Sab are determined by the S matrix bootstrap [31]; the only one we need 
later is that of the Ai — A2 scattering, which takes the form 



-] (-1 (-1 (- 

5 I J 15 I J 5 I J 15 



To have an unambiguous definition of the quantum numbers /j entering the Bethe-Yang equa- 
tions (2.14), it is convenient to define phase shift functions Sai, which are continuous and odd 
functions of tlie rapidity difference 0; we achieve this using tlie following convention: 



SabiB) = SabiO)e 



iSabid) 
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where Sab is uniquely specified by continuity and the branch choice 



SabiO) = 

and it is an odd function of 6 due to the following property of the scattering amplitude: 

Sabie)Sab{-9) = 1 

(which also implies Sab{0) = ±1)- The above redefinition of the phase shift compared to the 
original one in eqn. (2.15) contains as a special case the Lee- Yang definition (3.7) and also 
entails appropriate redefinition of quantum numbers depending on the sign of -S'o6(0). 

3.2.1 Truncated fermionic space approach for the Ising model 

The conformal Ising model can be represented as the theory of a massless Majorana fermion 
with the action ^ 

On the conformal plane the model has two sectors, with the mode expansions 



X^rez+i ^ Neveu-Schwarz (NS) sector 

^rez^rZ"'^"^^'^ Ramond (R) sector 



and similarly for the antiholomorphic field ip. The Hilbert space is the direct sum of a certain 
projection of the NS and R sectors, with the Virasoro content 

Hlsing = Vh^Vh 

(s) 

The spin field a connects the NS and R sectors, and its matrix elements B^j' in the sector with 
a given conformal spin s (of. eqn. (3.10)) can be most conveniently computed in the fermionic 
basis using the work of Yurov and Zamolodchikov [33], who called this method the truncated 
fermionic space approach. The fermionic basis can easily be chosen orthonormal, and thus 
in this case the metrics G^^^ on the spin subspaces (cf. eqn. (3.9)) are all given by unit 
matrices of appropriate dimension. Apart from the choice of basis all the calculation proceeds 
very similarly to the case of the Lee- Yang model. Energy and volume is measured in units of 
the lowest particle mass m = mi and using relation (3.14) one can write the dimensionless 
Hamiltonian in the form (3.8). The highest cutoff we use is Ccut = 30, in which case the 
Hilbert space contains around three thousand vectors (slightly depending on the value of the 
spin chosen). 

We remark that the energy density operator can be represented in the fermionic language 

as 

e = 'tptp 

which makes the evaluation of its matrix elements in the fermionic basis extremely simple. 
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3.2.2 Form factors of the energy density operator e 

The form factors of the operator e in the Eg model were first calculated in [34] and their 
determination was carried further in [35]. The exact vacuum expectation value of the field e 
is given by [36] 

{e)=eh\hf/^^ , €,, = 2.00314... 

or in terms of the mass scale m = nii 

(e) = 0.45475 ••• X m (3.17) 

The form factors are not known for the general n-particle case in a closed form, i.e. no 
formula similar to that in (3.12) exists. They can be evaluated by solving the appropriate 
polynomial recursion relations derived from the form-factor axioms. We do not present explicit 
formulae here; instead we refer to the above papers. For practical calculations we used the 
results computed by Delfino, Grinza and Mussardo, which can be downloaded from the Web 
in Mathematica format [37]. 

Our interest in the Ising model is motivated by the fact that this is the simplest model in 
which form factors of an operator different from the perturbing one are known, and also its 
spectrum and bootstrap structure is rather complex, both of which stands in contrast with 
the much simpler case of scaling Lee- Yang model. 



3.3 Evaluating matrix elements of a local operator O in TCSA 
3.3.1 Identification of multi-particle states 

Diagonalizing the TCSA Hamiltonian (3.8) yields a set of eigenvalues and eigenvectors at each 
value of the volume, but it is not immediately obvious how to select the same state at different 
values of the volume. Therefore in order to calculate form factors it is necessary to identify 
the states with the corresponding many-particle interpretation. 

Finding the vacuum state is rather simple since it is the lowest lying state in the spin-0 
sector and its energy is given by 

Eo{L) = BL + ... 

where the ellipsis indicate residual finite size effects decaying exponentially fast with volume 
L and B is the bulk energy density which in the models we consider is exactly known (3.5, 
3.15). One-particle states can be found using that their energies can be expressed as 

El'\L) = BL+^l^^y + ^" + --- 

again up to residual finite size effects where s is the spin of the sector considered and i is the 
species label (every sector contains a single one-particle state for each species). 

Higher multi-particle states can be identified by comparing the measured eigenvalues to the 
levels predicted by the Bethe-Yang equations. Fixing species labels ii, . . . ,z„ and momentum 
quantum numbers /i, . . . , /„, eqns. (2.14) can be solved to give the rapidities 6i, . . . ,9n of the 
particles as function of the dimensionless volume parameter I = mL. Then the energy of the 
multi-particle state in question is 



Elt^il"^ {L) = BL + Y, cosh ~ek + ... 

=1 
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which can be compared to the spectrum. 

For each state there exists a range of the volume, called the scaling region, where L is large 
enough so that the omitted residual finite size effects can be safely neglected and small enough 
so that the truncation errors are also negligible. More precisely, the scaling region for any 
quantity depending on the volume can be defined as the volume range in which the residual 
finite size corrections and the truncation errors are of the same order of magnitude; since both 
sources of error show a dependence on the state and the particular quantity considered (as 
well as on the value of the cutoff), so does the exact position of the scaling region itself. 

In the scaling region, we can use a comparison between the Bethe-Yang predictions and the 
numerical energy levels to sort the states and label them by multi-particle quantum numbers. 
An example is shown in figure 3.1, where we plot the first few states in the spin-0 sector of the 
scaling Lee- Yang model and their identification in terms of multi-particle states is given. In 
this case, the agreement with the predicted bulk energy density and the Bethe-Yang levels in 
the scaling region is better than one part in 10^ for every state shown (with the TCSA cutoff 
taken at Ccut = 30). 

3.3.2 Evaluation of matrix elements 

Suppose that we computed two Hamiltonian eigenvectors as functions of the volume L (labeled 
by their quantum numbers in the Bethe-Yang description (2.14), omitting the particle species 
labels for brevity): 

\{h,...,In})L = Y.'^iih,...,In;L)\i) 

i 

j 

in the sector with spin s and spin s', respectively. Let the inner products of these vectors with 
themselves be given by 

Af = J2^iih,...,ln;L)G\f^j{h,...,ln;L) 
W = ^*.(/(,...,4;L)G;f*,(l(,...,4;L) 

It is important that the components of the left eigenvector are not complex conjugated. In 
the Ising model we work in a basis where all matrix and vector components are naturally real. 
In the Lee- Yang model, the TCSA eigenvectors are chosen so that all of their components 
are either purely real or purely imaginary depending on whether the basis vector \i) is 
an element of the h = h = or the h = h = —1/5 component in the Hilbert space. It 
is well-known that the Lee- Yang model is non-unitary, which is reflected in the presence of 
complex structure constants as indicated in (3.2). This particular convention for the structure 
constants forces upon us the above inner product, because it is exactly the one under which 
TCSA eigenvectors corresponding to different eigenvalues are orthogonal. We remark that by 
redefining the structure constants and the conformal inner product it is also possible to use a 
manifestly real representation for the Lee- Yang TCSA (up to some truncation effects that lead 
to complex eigenvalues in the vicinity of level crossings [18]). Note that the above conventions 
mean that the phases of the eigenvectors are fixed up to a sign. 
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e(0 




|2,0,-2> 



one-particle 



Figure 3.1: The first 13 states in the finite volume spectrum of scaling Lee- Yang model. We 
plot the energy in units of m (with the bulk subtracted): e(Z) = {E{L) — BL)/m, against 
the dimensionless volume variable I = m,L. n-particle states are labeled by . . . ,In), where 
the Ik are the momentum quantum numbers. The state labeled |2, 1, —3) is actually two-fold 
degenerate because of the presence of | — 2, —1, 3) (up to a splitting which vanishes as e~', cf. 
the discussion in subsection 4.3). The dots are the TCSA results and the continuous lines are 
the predictions of the Bethe-Yang equations (2.14). The points not belonging to any of the 
Bethe-Yang lines drawn are two- and three-particle states which are only partly contained in 
the first 13 levels due to line crossings, whose presence is a consequence of the integrability of 
the model. 
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Let us consider a spinless primary field O with scaling weights Aq = A©, which can be 
described as the matrix 

0\^'''^ = {i\0{z,z)\j)\^^,^^ , \i) enTcs{s',e^^t) , li) G WTCs(^ecut) 

between the two truncated conformal space sectors. Then the matrix element of O can be 
computed as 

m-2^«({/(, . . . , /a|O(0, 0)|{/i, . . . , In})L = 

^) E ^.(^i' - ' lUL)0^t"'Ulu ...,In;L) (3.18) 

where the volume dependent prefactor comes from the transformation of the primary field O 
under the exponential map (3.1) and we wrote the equation in a dimensionless form using the 
mass scale m. The above procedure is a generalization of the one used by Guida and Magnoli 
to evaluate vacuum expectation values in [38]; it was extended to one-particle form factors in 
the context of the tricritical Ising model by Fioravanti et al. in [39]. 

4 Numerical results for elementary form factors 

4.1 Vacuum expectation values and one-particle form factors 
4.1.1 Scaling Lee- Yang model 

Before the one-particle form factor we discuss the vacuum expectation value. Let us define 
the dimensionless function 

0(Z) = -zm2/^(0|$|0)L 

where the finite volume expectation value is evaluated from TCSA using (3.18). We performed 
measurement of ij!) as a function of both the cutoff ecut = 21 ... 30 and the volume Z = 1 ... 30 
and then extrapolated the cutoff dependence fitting a function 

(where the exponent was chosen by verifying that it provides an optimal fit to the data). The 
data corresponding to odd and even values of the cutoff must be extrapolated separately [17], 
therefore one gets two estimates for the result, but they only differ by a very small amount (of 
order 10"^ at / = 30 and even less for smaller volumes). The theoretical prediction for (/>(Z) is 

(t){l) = 1.239394325 • ■ ■ + 0(6"') 

The numerical result (after extrapolation) is shown in figure 4.1 from which it is clear that 
there is a long scaling region. Estimating the infinite volume value from the flattest part of 
the extrapolated curve (at I around 12) we obtain the following measured value 

= oo) = 1.23938... 

where the numerical errors from TCSA are estimated to affect only the last displayed digit, 
which corresponds to an agreement within one part in 10^. 
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1.24 
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Figure 4.1: The vacuum expectation value of $ in finite volume. The dashed Une shows the 
exact infinite volume value, while the continuous line corresponds to eqn. (4.1). 

There is also a way to compute the leading exponential correction, which was derived by 
Delfino [40]: 

= ($) + 1 V F2{iTT,Q)uKo{mir) + ... (4.1) 

TT ^ — ' 

i 

where 

/•oo 

Ko{x) = de cosh0e-^^°"''^ 

is the modified Bessel-function, and the summation is over the particle species i (there is only 
a single term in the scaling Lee- Yang model). This agrees very well with the numerical data, as 
demonstrated in table 4.1 and also in figure 4.1. Using Liischer's finite- volume perturbation 
theory introduced in [21], the correction term can be interpreted as the sum of Feynman 
diagrams where there is exactly one propagator that winds around the cylinder, and therefore 
eqn. (4.1) can be represented graphically as shown in figure 4.2. 

To measure the one-particle form factor we use the correspondence (2.16) between the 
finite and infinite volume form factors to define the dimensionless function 

h{l) = -im^l' + {27:sff (0|$|M)l 

where |{s})l is the finite volume one-particle state with quantum number I = s i.e. from the 
spin-s sector. The theoretical prediction for this quantity is 

/f (0 = 1.0376434349 • • • + 0(6"') (4.2) 

The numerical results (after extrapolation in the cutoff) are shown in figure 4.3. The scaling 
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Figure 4.2: Graphical representation of eqn. (4.1). 
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(/>(/) (predicted) 


0(0 (TCSA) 


2 


1.048250 


1.112518 


3 


1.184515 


1.195345 


4 


1.222334 


1.224545 


5 


1.233867 


1.234396 


6 


1.237558 


1.237698 


7 


1.238774 


1.238811 


8 


1.239182 


1.239189 


9 


1.239321 


1.239317 


10 


1.239369 


1.239360 


11 


1.239385 


1.239373 


12 


1.239391 


1.239375 



Table 4.1: Comparison of eqn. (4.1) to TCSA data 
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Figure 4.3: One-particle form factor from sectors with spin s = 0,1,2. The continuous hne 
shows the exact infinite volume prediction. 



region gives the following estimates for the infinite volume limit: 

/0(/ = oo) = 1.037654... 
f 1(1 = oo) = 1.037650... 
fl[l = oo) = 1.037659... 

which show good agreement with eqn. (4.2) (the relative deviation is again around 10~^, as 
for the vacuum expectation value). 



4.1.2 Ising model in magnetic field 

For the Ising model, we again start with checking the dimensionless vacuum expectation value 
for which, using eqn. (3.17) we have the prediction 

m) = —{e)L = 0.45475 • • • + O f e"') 

where m = mi is the mass of the lightest particle and I = mL as before. The TCSA data 
are shown in figure 4.4. Note that there is substantial dependence on the cutoff Ccut and also 
that extrapolation in Ccut is really required to achieve good agreement with the infinite volume 
limit. Reading off the plateau value from the extrapolated data gives the estimate 

— (e) =0.4544... 

m 

for the infinite volume vacuum expectation value, which has 8 • 10~^ relative deviation from 
the exact result. Our first numerical comparison thus already tells us that we can expect much 
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Figure 4.4: Measuring the vacuum expectation value of e in the Ising model 



larger truncation errors than in the Lee- Yang case. It is also clear from figure 4.4 that in order 
to attain suitable precision in the Ising model extrapolation in the cutoff is very important. 
Defining the function 

m = {e)L/{e) 

we can calculate the leading exponential correction using eqn. (4.1) and the exact two-particle 
form factors from [37]. It only makes sense to include particles i = 1,2,3 since the contribution 
of the fourth particle is subleading with respect to two-particle terms from the lightest particle 
due to 7714 > 27771. The result is shown in figure 4.5; we do not give the data in numerical 
tables, but we mention that the relative deviation between the predicted and measured value 
is better than 10~^ in the range 5 < Z < 10. 

From now on we normalize all form factors of the operator c by the infinite volume vacuum 
expectation value (3.17), i.e. we consider form factors of the operator 

* = e/(e) (4.3) 

which conforms with the conventions used in [35, 37]. We define the dimensionless one-particle 
form factor functions as 

(2 \ 
+(27r.)M (0|^|{.}).L 

In the plots of figure 4.6 we show how these functions measured from TCSA compare to 
predictions from the exact form factors for particles i = 1, 2, 3 and spins s = 0, 1, 2, 3. 

It is evident that the scaling region sets in much later than for the Lee- Yang model; 
therefore for the Ising model we do not plot data for low values of the volume (all plots 
start from I ^ 10. ..15). This also means that truncation errors in the scaling region are also 
much larger than in the scaling Lee- Yang model; we generally found errors larger by an order 
of magnitude after extrapolation in the cutoff. We remark that extrapolation improves the 
precision by an order of magnitude compared to the raw data at the highest value of the cutoff. 

Note the rather large finite size correction in the case of A3. This can be explained 
rather simply as the presence of a so-called jU-term. We can again apply Liischer's finite- 
volume perturbation theory, which we use in the form given by Klassen and Melzer in [22] 
for finite volume mass corrections. The generalization to one-particle matrix elements is 



22 



i-O- 



1 

0.95 
0.9 

0.85 
0.8 

0.75 



~i 1 1 1 

^ * ^ 

f 

/ ■ extrapolated 

,' predicted 

/ infinite volume 



10 



Figure 4.5: The volume dependence of the vacuum expectation value of e in the Ising model, 
showing the extrapolated value and the prediction from eqn. (4.1), normalized by the value 
in the infinite volume limit. 

straightforward, and for a static particle it gives the diagram depicted in figure 4.7, whose 
contribution has the volume dependence 

Q-"^^^^ , = \Jml - ^ = 0.10453 • • • x mi 

i.e. we can expect a contribution suppressed only by e^*^'^'. A numerical fit of the /-dependence 
in the s = case is perfectly consistent with this expectation. As a result, no scaling region can 
be found, because truncation errors are too large in the volume range where the exponential 
correction is suitably small. We do not elaborate on this issue further here; we only mention 
that starting from this point there are other interesting observations that can be made, and 
we plan to return to them in a separate publication [26]. 

4.2 Two-particle form factors 
4.2.1 Scaling Lee- Yang model 

Following the ideas in the previous subsection, we can again define a dimensionless function 
for each two-particle state as follows: 

f2{l)hi, = -im^^H(i\^\{Ii,l2})L , l = mL 
Relation (2.16) gives the following prediction in terms of the exact form- factors: 

2/5 

h{l)i,i, = ''"^ Fi{h{l)M)) + 0(e-') (4.4) 

\l Piim),02{l)) 

where 0\{l)^ ^2(0 solve the Bethe-Yang equations 

/sinh^i + 5(^1 -6*2) = 27r/i 
Zsinh^2 + 5(^2 -6*1) = 27r/2 
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Figure 4.6: One-particle form factors measured from TCSA (dots) compared to the infinite- 
volume prediction from exact form factors. All numerical data have been extrapolated to 
Ccut = oo and s denotes the Lorentz spin of the state considered. The relative deviation in 
the scaling region is around 10~^ for Ai and A2; there is no scaling region for (see the 
discussion in the main text). 
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Figure 4.7: Leading finite size correction (a so-called //-term) to the one-particle form factor 

of ^3, which resiihs from the process of sphtting up into two copies of Ai which then wind 
around the cylinder once before recombining into again. 



and the density of states is given by 

Pii(^i,^2) = cosh 6*1 cosh 6*2 + / cosh 6'i(^(6'2 - ^i) +lcosh.02ip{0i - O2) 
the phase shift S is defined according to eqn. (3.7) and 

d5{e) 



dO 



There is a further issue to take into account: the relative phases of the multi-particle states are 
a matter of convention and the choice made in subsection 3.3.2 for the TCSA eigenvectors may 
differ from the convention adapted in the form factor bootstrap. Therefore in the numerical 
work we compare the absolute values of the functions /2(0 computed from TCSA with those 
predicted from the exact form factors. Note that this issue is present for any non-diagonal 
matrix element, and was in fact tacitly dealt with in the case of one-particle matrix elements 
treated in subsection 4.1.1. 

The prediction (4.4) for the finite vohime two-particle form factors is compared with spin-0 
states graphically in figure 4.8 and numerically in table 4.2, while the spin-1 and spin-2 case is 
presented in figure 4.9 and in table 4.3. These contain no more than a representative sample 
of our data: we evaluated similar matrix elements for a large number of two-particle states 
for values of the volume parameter I running from 1 to 30. The behaviour of the relative 
deviation is consistent with the presence of a correction of e~' type up to i ~ 9 ... 10 (i.e. the 
logarithm of the deviation is very close to being a linear function of I), and after Z ~ 16 ... 18 
it starts to increase due to truncation errors. This is demonstrated in figure 4.10 using the 
data presented in table 4.3 for spin-1 and spin-2 states^, but it is equally valid for all the other 
states we examined. In the intermediate region I ~ 10 ... 16 the two sources of numerical 
deviation are of the same order, and so that range can be considered as the optimal scaling 



^Note that the dependence of the logarithm of the deviation on the volume is not exactly linear because 
the residual finite size correction can also contain a factor of some power of I, and so it is expected that a logZ 



contribution is also present in the data plotted in figure 4.10. 
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Figure 4.8: Two-particle form factors in the spin-0 sector. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 

region: according to the data in the tables agreement there is typically around 10~^ (relative 
deviation). It is also apparent that scaling behaviour starts at quite low values of the volume 
(around / ~ 4 the relative deviation is already down to around 1%). 

It can be verified by explicit evaluation that in the scaling region the Betlie-Yang density of 
states (p) given in (2.17) differs by corrections of relative magnitude lO""-*^ — 10~^ (analytically: 
of order 1//) from the free density of states (p°) in (2.13), and therefore without using the 
proper interacting density of states it is impossible to obtain the precision agreement we 
demonstrated. In fact the observed 10"'^ relative deviation corresponds to corrections of order 

at I = 10, but it is of the order of estimated truncation errors^. 

These results are very strong evidence for the main statement in (4.4) (and thus also 
(2.16)), namely, that all l/L corrections are accounted by the proper interacting state density 
factor and that all further finite size corrections are just residual finite size effects decaying 
exponentially in L. In section 4.3 we show that data from higher multi-particle form factors 
fully support the above conclusions drawn from the two-particle form factors. 

4.2.2 Ising model in magnetic field 

In this case, there is some further subtlety to be solved before proceeding to the numerical 
comparison. Namely, there are spin-0 states which are parity refiections of each other, but 
are degenerate according to the Bethe-Yang equations. An example is the state |{1,— l})i2 
in figure 4.11 (b), which is degenerate with |{ — l,l})i2 to all orders in l/L. In general the 
degeneracy of these states is lifted by residual finite size effects (more precisely by quantum 
mechanical tunneling - a detailed discussion of this mechanism was given in the framework of 

^Truncation errors can be estimated by examining the dependence of the extracted data on the cutoff ecut, 
as well as by comparing TCSA energy levels to the Bethe-Yang predictions. 
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h = 1/2 , h = -1/2 


h = 3/2 , h = -3/2 


h = 5/2 , h = -5/2 
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TCSA 
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TCSA 


FF 


TCSA 


FF 


2 


0.102780 


0.120117 


0.058158 


0.066173 


0.039816 


0.045118 


4 


0.085174 


0.086763 


0.058468 


0.059355 


0.042072 


0.042729 
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0.056828 


0.056769 


0.050750 


0.050805 


0.039349 


0.039419 
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0.036058 


0.035985 


0.042123 


0.042117 


0.035608 


0.035614 


10 


0.023168 


0.023146 


0.034252 


0.034248 


0.031665 


0.031664 


12 


0.015468 


0.015463 


0.027606 


0.027604 


0.027830 


0.027828 


14 


0.010801 


0.010800 


0.022228 


0.022225 


0.024271 


0.024267 


16 


0.007869 


0.007867 


0.017976 


0.017972 


0.021074 


0.021068 


18 


0.005950 


0.005945 


0.014652 


0.014645 


0.018268 


0.018258 


20 


0.004643 


0.004634 


0.012061 


0.012050 


0.015844 


0.015827 



Table 4.2: Two-particle form factors |/2(0I ™ the spin-0 sector 
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0.054711 


0.064623 


0.074763 


0.042031 


0.047672 


4 


0.072104 


0.073571 


0.049790 


0.050566 


0.062533 


0.063932 


0.044034 


0.044716 


6 


0.056316 


0.056444 


0.045031 


0.045100 


0.051828 


0.052009 


0.040659 


0.040724 


8 


0.042051 


0.042054 


0.039191 


0.039193 


0.041370 


0.041394 


0.036284 


0.036287 


10 


0.031146 


0.031144 


0.033469 


0.033467 


0.032757 


0.032759 


0.031850 


0.031849 


12 


0.023247 


0.023245 


0.028281 


0.028279 


0.026005 


0.026004 


0.027687 


0.027684 


14 


0.017619 


0.017616 


0.023780 


0.023777 


0.020802 


0.020799 


0.023941 


0.023936 


16 


0.013604 


0.013599 


0.019982 


0.019977 


0.016808 


0.016802 


0.020659 


0.020652 


18 


0.010717 


0.010702 


0.016831 


0.016822 


0.013735 


0.013724 


0.017835 


0.017824 


20 


0.008658 


0.008580 


0.014249 


0.014227 


0.011357 


0.011337 


0.015432 


0.015413 



Table 4.3: Two-particle form factors |/2(0I the spin-1 and spin-2 sectors 
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Figure 4.10: Estimating the error term in (4.4) using the data in table 4.3. The various plot 
symbols correspond to the same states as specified in figure 4.9. 
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the fc-folded sine-Gordon model in [41]). Since the finite volume spectrum is parity symmetric, 
the TCSA eigenvectors correspond to the states 

|{1, -1})± = i= (|{1, -1})12,L ± |{-1, 1})i2,l) 

Because the Hilbert space inner product is positive definite, the TCSA eigenvectors |{1, — l})f2 
can be chosen orthonormal and the problem can be resolved by calculating the form factor 
matrix element using the two-particle state vectors 



^(i{i>-i})r2,L±ia-i})r2,L) 



V2 

Because the Ising spectrum is much more complicated than that of the scaling Lee- Yang model 
(and truncation errors are larger as well), we only identified two-particle states containing two 
copies of Ai , or an Ai and an A2. The numerical results are plotted in figures 4.11 (a) and 
(b), respectively. The finite volume form factor functions of the operator ^ (4.3) are defined 
as 



where 



hi (0/,/, = V/^ii(^i(0>^2(0)(o|^'|{/i,/2})ii 



/ sinh^i + 5u{9i - 62) = l-Kh 
lsmhe2 + Sn{02-0i) = 2Trl2 

Pii(^i5^2) = cosh ^1 cosh ^2 + Zcosh^i (^11(^2 — ^1) + / cosh02</'ii(^i — ^2) 

dSn{9) 



dd 



and 



with 



/12 {i)i,i, = \^piiiOi,02){omh,h})i2 



Zsinheii + SuiOi - 62) = 27r/i 

— Zsinh^2 + '^12(^2 - ^1) = 27r/2 
mi 

Pi2{0i,92) = —/^ cosh ^1 cosh ^2 + ^ cosh ^1(^12 (6*2 - ^i) + cosh 6'2V'i2 (6*1 - O2) 
mi mi 

and are compared against the form factor functions 

F^{eiil),92{l))u 

and 

F^{ei{i),e2{i))i2 

respectively. 

Although (as we already noted) truncation errors in the Ising model are much larger than 
in the Lee- Yang case, extrapolation in the cutoff improves them by an order of magnitude 
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compared to the evaluation at the highest cutoff (in our case 30). After extrapolation, devia- 
tions in the scaling region become less than 1% (with a minimum of around 10~^ in the ^i^i, 
and 10~^ in the A1A2 case), and even better for states with nonzero total spin. As noted 
in the previous subsection this means that the numerics is really sensitive to the dependence 
of the particle rapidities and state density factors on the interaction between the particles; 
generally the truncation errors in the extrapolated data are about two orders of magnitude 
smaller than the interaction corrections. 

It is a general tendency that the agreement is better in the sectors with nonzero spin, and 
the scaling region starts at smaller values of the volume. This is easy to understand for the 
energy levels, since for low-lying states nonzero spin generally means higher particle momenta. 
The higher the momenta of the particles, the more the Bethe-Yang contributions dominate 
over the residual finite size effects. This is consistent with the results of Rummukainen and 
Gottlieb in [42] where it was found that resonance phase shifts can be more readily extracted 
from sectors with nonzero momentum; our data show that this observation carries over to 
general matrix elements as well. 

4.3 Many-particle form factors 
4.3.1 Scaling Lee- Yang model 

We also performed numerical evaluation of three and four-particle form factors in the scaling 
Lee- Yang model; some of the results are presented in figures 4.12 and 4.13, respectively. For 

the sake of brevity we refrain from presenting explicit numerical tables; we only mention 
that the agreement between the numerical TCSA data and the prediction from the exact form 
factor solution is always better than 10~^ in the scaling region. For better visibility we plotted 
the functions 

fk{l)h...h = -im^^''\/pk{Oi, ek)m\{h, . . .,h})L , l = mL 

for which relation (2.16) gives: 

h{l)h...iu=-im^'''Ft{eu...A) + 0{e-') (4.5) 

Due to the fact that in the Lee- Yang model there is only a single particle species, we introduced 
the simplified notation pn for the n-particle Jacobi determinant. 

The complication noted in subsection 4.2.2 for the Ising state |{1,— l})i2 is present in 
the Lee- Yang model as well. The Bethe-Yang equations give degenerate energy vahies for 
the states . . . ,Ik}) l and |{— /fc, . . . , — (as noted before, the degeneracy is lifted by 
quantum mechanical tunneling). For states with nonzero spin this causes no problem, because 
these two states are in sectors of different spin (their spins differ by a sign) and similarly there 
is no difficulty when the two quantum number sets are identical, i.e. 

{^1, ■■■ .Ik} = {—h, ■■■ , —Ik} 

since then there is a single state. However, there are states in the zero spin sector (i.e. with 
Y^f^ Ik = 0) for which 

{Ii,...,Ik}^{-Ii,---,-Ik} 

We use two such pairs of states in our data here: the three-particle states 1{3, — 1, — 2})^, 
|{2, 1, -3}) L and the four-particle states |{7/2, 1/2, -3/2, -5/2}) l , |{5/2, 3/2, -1/2, -7/2}) l ■ 
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Figure 4.11: Two-particle form factors in the Ising model. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 
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Again, the members of such pairs are related to each other by spatial reflection, which is a sym- 
metry of the exact finite- volume Hamiltonian and therefore (supposing that the eigenvectors 
are orthonormal) the finite volume eigenstates correspond to 

. . . , Ik})t = . . . , Ik})L ± \{-Ik, . . . , -h})L) 

and this must be taken into account when evaluating the form factor matrix elements. In 
the Lee- Yang case, however, the inner product is not positive definite (and some nonzero 
vectors may have zero "length", although this does not happen for TCSA eigenvectors, because 
they are orthogonal to each other and the inner product is non-degenerate), but there is a 
simple procedure that can be used in the general case. Suppose the two TCSA eigenvectors 
corresponding to such a pair are vi and V2- Then we can define their inner product matrix as 

using the TCSA inner product (3.9). The appropriate basis vectors of this two-dimensional 
subspace, which can be identified with . . . , Ik})L and \{—Ik, ■ ■ ■ , can be found by 

solving the two-dimensional generalized eigenvalue problem 

g ■ w = XP ■ w 

for the vector (^1,1^2) describing orientation in the subspace, with 



P = 



1 

1 



This procedure has the effect of rotating from the basis of parity eigenvectors to basis vectors 
which are taken into each other by spatial reflection. 



4.3.2 Ising model in a magnetic field 

As we already noted, it is much harder to identify^ higher states in the Ising model due to the 
complexity of the spectrum, and so we only performed an analysis of states containing three 
Ai particles. We define 

/ill = ^Jpium),e2ii),e3ii))m\{h,i2,i3})ui 

where 6i{l) are the solutions of the three-particle Bethe-Yang equations in (dimensionless) 
volume I and pm is the appropriate 3-particle determinant. The results of the comparison 
can be seen in figure 4.14. The numerical precision indicated for two-particle form factors at 
the end of subsection 4.2.2, as well as the remarks made there on the spin dependence apply 
here as well; we only wish to emphasize that for AiAiAi states with nonzero total spin the 
agreement between the extrapolated TCSA data and the form factor prediction in the optimal 
part of the scaling region is within 2 x 10~^. 

■*To identify ^i^i^i states it is necessary to use at least Ccut = 22 or 24 and even then the agreement with 
the Bethe-Yang prediction is still only within 20%, but the identification can be made for the first few AiAiAi 
states using data up to Ccut = 30. Truncation errors are substantially decreased by extrapolation to Ccut = 00. 
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Figure 4.12: Three-particle form factors in the spin-0 sector. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 
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Figure 4.13: Four-particle form factors in the spin-0 sector. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 
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Figure 4.14: Three-particle form factors in the Ising model. Dots correspond to TCSA data, 
while the lines show the corresponding form factor prediction. 



5 General form factors without disconnected pieces 

Let us consider a matrix element of the form 

n...jm{{l'i, ■ ■ ■ , l'm}\O{0, 0)|{/i, . . . , /n})n...«„,L 

Disconnected pieces are known to appear when there is at least one particle in the state on 
the left which occurs in the state on the right with exactly the same rapidity. The rapidities 
of particles as a function of the volume are determined by the Bethe-Yang equations (2.14) 

QkjOi, ■ ■ ■,On)ii...i„ = rwjfeLsinhgfc + y^(^ifeii(4 - Oi) = 27rlfc , k = l,...,n 

l^k 

and 

Qkie[, O'Jn-jru = m,,L sinh 9', + ^ dj.j, (6', - O'l) = 27r4 , k = l,...,m 

Due to the presence of the scattering terms containing the phase shift functions 6, equality of 
two quantum numbers and I'l does not mean that the two rapidities themselves are equal 
in finite volume L. It is easy to see that there are only two cases when exact equality of some 
rapidities can occur: 

1. The two states are identical, i.e. n = m and 

{il • • • jm} = {k--- in} 
{-^1) • • • = {h,---,In} 

in which case all the rapidities are pairwise equal, or 
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2. Both states are parity symmetric states in the spin zero sector, i.e. 

{/!,...,/„} = {-/„,..., -/i} 
{-^1 ) • • • ) -^m} — { -^mi • • • 1 -^l} 

and the particle species labels are also compatible with the symmetry, i.e. in+i-k = ik 
and jm+i-k = jk- Furthermore, both states must contain one (or possibly more, in a 
theory with more than one species) particle of quantum number 0, whose rapidity is then 
exactly for any value of the volume L due to the symmetric assignment of quantum 
numbers. 

Discussion of such matrix elements raises many interesting theoretical considerations and is 
postponed to the followup paper [20]; here we only concentrate on matrix elements for which 
there are no disconnected contributions. 

5.1 Scaling Lee- Yang model 

In this model there is a single particle species, so we can introduce the following notations: 
fknilfci = ■ ■ . ,411^(0, 0)|{/i, . . . 

and also 

fkn{l)t::,l = -im'/'^Pk{0[, K)^JpniOu en){{l[, • • • , 411^(0, 0)|{/l, . . . , In}) L 

for which relation (2.20) yields 

^Jpniel,...,en)Pk{0[,■■■,O'm) 

fkn{iyt::fn = -im''/'F^+n0'k + iT^,---,O'i+i^,0i,---,On) + O{e-') (5.1) 

For the plots we chose to display / or / depending on which one gives a better visual picture. 
The numerical results shown here are just a fraction of the ones we actually obtained, but 
all of them show an agreement with precision 10~^ — lO^"* in the scaling region (the volume 
range corresponding to the scaling region typically varies depending on the matrix element 
considered due to variation in the residual finite size corrections and truncation effects). 

The simplest cases involve one and two-particle states: the one-particle-one-particle data 
in figure 5.1 actually test the two-particle form factor F.f , while the one-particle-two-particle 
plot 5.2 corresponds to F.f (we obtained similar results on using matrix elements /22)- 
Note that in contrast to the comparisons performed in subsections 4.2 and 4.3, these cases 
involve the form factor solutions (3.12) at complex values of the rapidities. In general, all 
tests performed with TCSA can test form factors at rapidity arguments with imaginary parts 
or TT, which are the only parts of the complex rapidity plane where form factors eventually 
correspond to physical matrix elements. 

One-particle-three-particle and one-particle-four-particle matrix elements /13 and /14 con- 
tribute another piece of useful information. We recall that there are pairs of parity-related 
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Figure 5.1: One-particle-one-particle form factors in Lee- Yang model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 5.2: One-particle-two-particle form factors in Lee- Yang model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 
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({1}|$|{1,0-1}) 
({l}|$|{2,0-2}) 
({0}|$|{3-l-2}> 
{{l}|$|{3-l-2}) 

{{1}|$|{-3,1,2}) 



Figure 5.3: One-particle-three-particle form factors in Lee- Yang model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 

states in the spin-0 factors which we cannot distinguish in terms of their elementary form 
factors. In subsection 4.3 we showed the example of the three-particle states 

|{3,-l,-2})i and |{2,1,-3})l 

and the four-particle states 

|{7/2,l/2,-3/2,-5/2})i and |{5/2, 3/2, -1/2, -7/2})^ 

In fact it is only true that they cannot be distinguished if the left state is parity-invariant. 
However, using a one-particle state of nonzero spin on the left it is possible to distinguish and 
appropriately label the two states, as shown in figures 5.3 and 5.4. This can also be done using 
matrix elements with two-particle states of nonzero spin: the two-particle-three-particle case 
/23 is shown in 5.5 (similar results were obtained for 724). Examining the data in detail shows 
that the identifications provided using different states on the left are all consistent with each 
other. 

It is also interesting to note that the /14 (figure 5.4) and /23 data (figure 5.5) provide a test 
for the five-particle form factor solutions F5. This is important since it is progressively harder 
to identify many-particle states in the TCSA spectrum for two reasons. First, the spectrum 
itself becomes more and more dense as we look for higher levels; second, the truncation errors 
grow as well. Both of these make the identification of the energy levels by comparison with 
the predictions of the Bethe-Yang equations more difficult; in the Lee- Yang case we stopped 
at four-particle levels. However, using general matrix elements and the relations (5.1) we can 
even get data for form factors up to 8 particles, a sample of which is shown in figures 5.6 (/33 
and /44, corresponding to 6 and 8 particle form- factors) and 5.7 (734 which corresponds to 7 
particle form factors). 

5.2 Ising model in magnetic field 

In the case of the Ising model, we define the functions 

fji-jm-M-iM = \l Pii-in{dl, ■ ■ ^n)Pji...jm{^'i, • • • , 6'to)x ji...jm(Ul' ■ ■ ■ >-^m}l*IUl' ■ ■ ■ > -fn})n...i„,L 
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Figure 5.4: One-particle-four-particle form factors in Lee- Yang model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 5.5: Two-particle-three-particle form factors in Lee- Yang model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 5.6: Three-particle-three-particle and four-particle-four-particle form factors in Lee- 
Yang model. Dots correspond to TCSA data, while the lines show the corresponding form 
factor prediction. 
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Figure 5.7: Three-particle-four-particle form factors in Lee- Yang model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 



39 



which are compared against form factors 

^m+n0m + IT^ , ■ ■ ■ , O'l + ITT , 9i, . . . , On)jm...jiii...i„ 

where 6i and 9'^ denote the rapidities obtained as solutions of the appropriate Bethe-Yang 
equations at the given value of the volume. We chose states for which the necessary form 
factor solution was already known (and given in [37]) i.e. we did not construct new form 
factor solutions ourselves. 

One-particle-one-particle form factors are shown in figure 5.8; these provide another nu- 
merical test for the two-particle form factors examined previously in subsection 4.2.2. One- 
particle-two-particle form factors, besides testing again the three-particle form factor yli^i^i 
(figure 5.9 (a)) also provide information on A1A1A2 (figure 5.9 (b)). 

Finally, one-particle-three-particle and two-particle-two-particle matrix elements can be 
compared to the ^lAi^diTli form factor, which again shows that by considering general matrix 
elements we can go substantially higher in the form factor tree than using only elementary 
form factors. 

We remark that the cusps on the horizontal axis in the form factor plots correspond to 
zeros where the form factors change sign; they are artifacts introduced by taking the absolute 
value of the matrix elements. The pattern of numerical deviations between TCSA data and 
exact form factor predictions is fully consistent with the discussion in the closing paragraphs 
of subsections 4.2.2 and 4.3.2. The deviations in the scaling region are around 1% on average, 
with agreement of the order of 10~^ in the optimal range. 

6 Conclusions 

In this work we gave an expression for the matrix elements of local operators between general 
multi-particle states in finite volume which is valid to any order in the expansion in the inverse 
of the volume 1/L. It was shown in section 2.2 using a very general argument that all the 
remaining volume dependence is non-analytic in 1/L (it is given by residual finite size effects 
vanishing exponentially with increasing volume). It is also clear that the derivation itself 
does not depend on integr ability, neither it is restricted to 1 -|- 1 dimensional field theories 
and therefore relations (2.16) and (2.20) can be extended to general quantum field theories 
(substituting the rapidities with appropriate kinematical parametrization and the pn with the 
proper state densities), with the only condition that their spectrum of excitations must possess 
a mass gap. 

1-i-l dimensional integrable field theories are special in the respect that multi-particle states 
in finite volume can be described using the Bethe-Yang equations (2.14) and so the n-particle 
state density Pn can be obtained in the general closed form (2.17). Another important feature is 
that there are exact results for matrix elements of local operators in infinite volume which can 
be obtained from the form factor bootstrap briefly reviewed in section 2.1. Therefore they are 
ideal toy models to test ideas about finite size corrections. Such an approach is also interesting 
due to a fundamental property of the bootstrap, namely that it is only indirectly related to 
the actual Lagrangian (or Hamiltonian) field theory. As we discussed in the introduction, 
testing the conjectured form factors against field theory usually involves calculating two-point 
functions using spectral representations, or sum rules derived from such expansions; however, 
direct non-perturbative comparison of the actual form factors to matrix elements computed 
from the field theory have been very restricted so far. 
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(a) Ai — Ai matrix elements 
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(b) Ai — A2 matrix elements 

Figure 5.8: One-particle-one-particle form factors in the Ising model. Dots correspond to 
TCSA data, while the lines show the corresponding form factor prediction. 
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(a) A\ — A\A\ matrix elements 
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(b) A\ — A1A2 matrix elements 

Figure 5.9: One-particle-two-particle form factors in the Ising model. Dots correspond 
TCSA data, while the lines show the corresponding form factor prediction. 
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Figure 5.10: One-particle-three-particle and two-particle-two-particle form factors in the Ising 
model. Dots correspond to TCSA data, while the lines show the corresponding form factor 
prediction. 
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Using TCSA we were able to give an extensive and direct numerical comparison between 
bootstrap results for form factors and matrix elements evaluated non-perturbatively. One of 
the advantages is that we can compare matrix elements directly, without using any proxy (such 
as a two-point function or a sum rule); the other is the very high precision of the comparison 
and also that it is possible to test form factors of many particles which have never been 
tested using spectral sums, mostly due to the fact that usually their contribution to spectral 
expansions is extremely small, and evaluating it also involves calculating multidimensional 
integrals to very high precision, which a numerically difhcult task. The second problem is 
actually related to the first, since due to the smallness of the contribution from higher particle 
terms all the lower ones must be evaluated to sufficiently high precision. Our approach, in 
contrast, makes it possible to have a test of entire one-dimensional sections of the form factor 
functions using the volume as a parameter, and the number of available sections only depends 
on our ability to identify multi-particle states in finite volume. 

Our results can also be viewed in the context of finite volume form factors [15] (which is 
also related to the problem of finite temperature form factors; for a review on the latter see [16] 
and references therein). The relations (2.16, 2.20) give finite volume form factors expressed 
with their infinite volume counterparts to all orders in 1/L (where L denotes the volume), i.e. 
up to exponentially decaying terms in L. This gives finite volume form factors in large volume 
with very high precision. On the other hand, what we determine numerically in TCSA are 
actually the finite volume form factors themselves, which is an approach that primarily works 
in small enough volume due to the truncation errors. In the Lee- Yang case, the combination 
of the two approaches gives the finite volume form factors involving up to four particles with 
better than 10~^ relative precision, as demonstrated by the excellent agreement in the scaling 
region of TCSA where their domains of validity overlap. For the Ising model the numerical 
precision is not as good, but with some care a precision of around 10~^ can be achieved for 
most of the matrix elements considered in this paper. 

An open question which is not discussed in this paper is the case of matrix elements with 
disconnected pieces. Results on such matrix elements are already available, but we postpone 
them to a followup paper [20] where we also plan to discuss many theoretical issues related to 
crossing and disconnected contributions in finite volume. 
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